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Abstract 

Using the transfer matrix method and the effective-mass approximation, the effect of resonant 
states on spin transport is studied in ZnSe/ZnMnSe/ZnSe/ZnMnSe/ZnSe structures under the 
influence of both electric and magnetic fields. The numerical results show that the ZnMnSe layers, 
which act as spin filters, polarize the electric currents. Variation of thickness of the central ZnSe 
layer shifts the resonant levels and exhibits an oscillatory behavior in spin current densities. It 
is also shown that the spin polarization of the tunneling current in geometrical asymmetry of the 
heterostructure where two ZnMnSe layers have different Mn concentrations, depends strongly on 
the thickness and the applied bias. 
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I. INTRODUCTION 



Tunneling and resonant tunneling processes involving electron spin show prominent per- 
spectives for new field spin-dependent fast electronics. Besides important applications in 
new spin-based multifunctional devices such as spin-field-effect transistor and spin-light- 
emitting diode, the spin- dependent resonant tunneling (SDRT) effect can also help us to 
more deeply understand the role of spin degree of freedom of the tunneling electron and the 
quantum size effects on spin transport processes m recent years, some attempts have 
been made to study the effects of SDRT in magnetic tunnel junctions in which the electrons 
tunnel from one ferromagnetic metal electrode to the other through a nonmagnetic metal 

layer and a thin insulator layer DQDDD 

On the other hands, spin-injection into semiconductors by using magnetic semiconductors 
has also been demonstrated, openingthe way to new all-semiconductor devices. II- VI diluted 
magnetic semiconductors (DMSs) |7j are known to be good candidates for effective spin 
injection into a nonmagnetic semiconductor (NMS) because their spin polarization is nearly 
100% and their conductivity is comparable to that of typical NMS. Moreover, II- VI DMSs 
can be n-type doped, thus avoiding the very fast spin precession that limits the applicability 
of III-V DMSs as spin injectors. A very promising II- VI DMS for spin injection is (Zn,Mn)Se, 
which has been previously used for spin injection experiments into GaAs 0, and ZnSe |^|. 
More recently, Slobodskyy et al. nj] using II- VI semiconductor layers, fabricated a magnetic 
resonant tunneling diode. Their magnetic device is based on a quantum well made of diluted 
magnetic semiconductor ZnMnSe between two ZnBeSe barriers and surrounded by highly 
n-type ZnSe layers. 

12 1 and Chang et al. 13] investigated theoretically the spin- 



Recently, Guo et al. 
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dependent transport in ZnSe/Zn 1 „ x Mn x Se DMS double barrier structures. The results 
showed that the spin polarization of the tunneling electrons can be tuned by changing 
the external magnetic and electric fields. The effect of SDRT on spin transport has also 



been studied in spin filter tunnel junctions I i|, 3,0, 12, G] in which a ferromagnetic 



semiconductor (FMS) such as EuS is used as a magnetic barrier. The large spin polarization 



achievable using magnetic barriers |19|, makes spin filtering a nearly ideal method for spin 
injection into semiconductors, enabling novel spintronic devices j^j. 



In Ref. [la, using two FMS layers, we examined the SDRT in a double barrier junction 
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and in the coherent tunneling regime. In the present work, we study theoretically the effects 
of SDRT in ZnSe/ZnMnSe/ZnSe/ZnMnSe/ZnSe structures. The left and right ZnSe layers 
are considered as emitter and collector attached to external leads. We assume that the car- 
rier wave vector parallel to the interfaces and the carrier spin are conserved in the tunneling 
process through the whole system. The first assumption is well justified for interfaces be- 
tween materials whose lattice constants are nearly equal; and the second one is also justified 
in our structure, because the sample dimensions are much smaller than the spin coherence 
length. Using the transfer matrix method and the nearly-free-electron approximation we will 
study the effects of resonant tunneling on spin- dependent current densities and the degree 
of spin polarization. 

The paper is organized as follows. In section II, the model is described and then the 
spin current densities and spin polarization for ZnSe/ZnMnSe heterostructures with two 
paramagnetic layers, are formulated. In section III, the numerical results for the symmetric 
and asymmetric structures are discussed in terms of the thickness of the ZnSe layer, which has 
been sandwiched between two paramagnetic layers. The results of this work are summarized 
in section IV. 

II. MODEL AND FORMALISM 

We consider a spin unpolarized electron current injected into 
ZnSe/Zn 1 _ :E Mn :E Se/ZnSe/Zn 1 _j / Mn ?/ Se/ZnSe structures in the presence of magnetic 
and electric fields along the growth direction (taken as z axis). The conduction electrons 
that contribute to the electric currents, interact with the 3d 5 electrons of the Mn ions via 
the sp-d exchange interaction. Due to the sp-d exchange interaction, the external magnetic 
field gives rise to the spin splitting of the conduction band states in the paramagnetic 
layers. Therefore, the injected electrons see a spin-dependent potential. In the framework 
of the parabolic-band effective mass approximation, the one-electron Hamiltonian of such 
system can be written as 

1 pV 7 

H = ^(P + eA) 2 + V S + V (z) + K z (z) - , (1) 

where the electron effective mass m* is assumed to be identical in all the layers, and the 
vector potential is taken as A = (0,Bx, 0). Here, V s = \g s HB<J • B describes the Zeeman 
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splitting of the conduction electrons, where a is the conventional Pauli spin operator. Vq(z) 
is the conduction band offset under zero magnetic field, which is the difference between the 
conduction band edge of the ZnMnSe layers and that of the ZnSe layer; V Uz (z) is the sp-d 
exchange interaction between the injected electron and the Mn ions and can be calculated 
within the mean field approximation. Hence, the sum of the two terms can be written as 



where, V 0L (r), N a L{R ) and (S zL (r)) = SB s [5ii B B/k B (T + T 0L{R) )] are respectively the con- 
duction band offset, the electronic sp-d exchange constant and the thermal average of zth 
component of Mn 2+ spin in the left (right) DMS layer. Here, Bs{x) is the Brillouin function 
and S = 5/2 is the spin of the Mn ions. a z = ±1/2 (or f, j) are the electron spin com- 
ponents along the magnetic field, x e ff = x(l — x) 12 and y e ff = y(l — y) 12 are the effective 
Mn concentrations due to the antiferromagnetic Mn-Mn coupling, while x and y are real 
Mn concentrations. Q(z) is the step function, L\ and L3 are respectively the widths of 
left and right DMS layers, and L2 is the width of the middle ZnSe layer. The last term 
in Eq. (JTJ) denotes the effect of an applied bias V a along the z axis on the system, where 
L = L\ + L 2 + L 3 is the total length of the considered structure along the growth direction. 
It is important to note that, our sample dimensions are much smaller than the spin coherent 
length in the semiconductors. Therefore we have neglected the effects of spin-flip processes 
in the Hamiltonian of the system. 

Here, we would like to point out that, due to the absence of any kind of scattering center 
for the electrons, the motion along the z-axis is decoupled from that of the x-y plane which 
is quantized in the Landau levels with energies E n = {n + l/2)fiu) C) where n = 0, 1, 2, • • • 
and uj c = eB/m*. In such case, the motion of electrons along the z axis can be reduced to 
a one- dimensional Schrodinger equation 



V (z) + V a .(z) = [V 0L -N a L a z x eff (S zL )]e(z)e(L 1 -z) 

+ [V 0R - N a R cr z y eff (S zR )]Q(z -L x - LjQfa + L 2 + L 3 - z) , (2) 



h 2 d 2 ^ az (z) 



+ V s + V (z) + V ff „(z) 



eV a z~ 



(3) 



2m* dz 2 



L 



where E z is the longitudinal energy of electrons traversing the heterostructure. 
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The general solution to the above Schrodinger equation is as follows: 



A laz e ik ^ z + B laz e- ikl "* z , z<0, 

A 2az Ai[p az (z)} + B 2(7z Bi[p az (z)}, 0<z<L x , 

A 3az Ai[p az {z)} + B 3az Bi[p az {z)], L 1 <z<L 1 + L 2 , 

A 4az Ai[p az (z)} + B Auz Bi[ P(7z {z)}, Li + L 2 < z < L, 



(4) 



.4 



ik h „,z 



+ B 5rT e- ik ^ z , 



z > L, 



where the coefficients Aj az and Bj az (with j—1-5) are constants which can be determined 
by the boundary conditions, A\[p az (z)] and Bi[p CTz (z)] are Airy functions, and 



k laz = J2m*(E z - V s )/h 



P*X Z ) 



k 5az = j2m*(E z -V s + eV a )/h 



V s + V (z) + V az (z)-E z - 



eV a X 



eV a z 



with 



(5) 

(6) 
(7) 

2m eV a ' (M 
The wave functions and their first derivatives in the five regions are matched at the 
interfaces between the regions. The matching results in a system of equations, which can be 
represented in a matrix form j^jj, 

/ 



A 



-h 2 L 



1/3 




total 



A 



OCT- 



(9) 



where M tota i is the transfer matrix that connects the incidence and transmission amplitudes 
and has the following form 



M total = M^(0)M 2 (0)M 2 \L 1 )M 3 (L 1 )M 3 -\L 2 )M 4 (L 2 )M^ 1 (L 3 )M 5 (L 3 ) . (10) 



Here, 



Mj(zi) 



(11) 



dz dz 

where, ipj{z) and ipj{z) are respectively the first and second term of the wave functions in 
each layer, without considering their coefficients. Since there is no reflection in region 5, the 
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coefficient B^ az in Eq. (jlj) is zero and the transmission coefficient of the spin a z electron, 
which is defined as the ratio of the transmitted flux to the incident flux, can be written as 

2 



T a ,(E z ,B,V a ] 



fa<7 z 



M tota i 



(12) 



where M^ tal is the left-upper element of the matrix M tota i defined in Eq. (jTDjl. It should 
be noted that the transmission coefficients depend on the incident energy E z , the magnetic 
field B, the applied bias V a and the spin orientation. The spin-dependent current density is 
connected with the transmission coefficients via 



00 r+00 

J„ z (B,V a ) = JoBY, T az (E z ,B,V a ) 



x{f[E z + (n + ^)hu + V a ] - f[E z + (n + ^)Hoj + V S + eV a }}dE z , (13) 

where Jo = e 2 /4ir 2 h 2 and f(E) = 1/ {l+exp[(E — Ep) / kgT]} is the Fermi-Dirac distribution 
function in which Ep denotes the emitter Fermi energy. 

The degree of spin polarization for electrons traversing the heterostructure is defined by 

n Ji(B,V a )-MB,V a ) 



j ; (s,y ) + j T (s,K) 

where Jj (J|) is the spin-up (spin-down) current density. 



(14) 



III. RESULTS AND DISCUSSION 

From Eq. (|13|) . we can evaluate numerically the spin-dependent current densities as a 
function of the thickness of the central ZnSe layer for the symmetric and the asymmetric 
structures, depending on the Mn concentration in both paramagnetic layers. In all of the 
numerical calculations we use m* = 0.16 m e (m e is the mass of the free electron), T = 
4.2 K, B = 4 T, Li = L 3 = 75 nm, g s = 1.1, and Ep = 5 meV. Figs. 1(a) and 1(b) 
show the thickness dependence of the current densities for spin-up and spin-down electrons 
respectively, in a symmetric (x = y = 0.05) structure. In both figures, Nqq:l = Noan = 0.26 
eV, Vol = Vqr = meV, T 0L = T 0R = 1.7 K [22]. By applying an external magnetic 
field, the sp-d exchange interaction gives rise to a giant spin splitting which is much larger 
than the Zeeman splitting of the conduction electrons. In such case, the degeneracy of 
the spin-up and spin-down electron states is lifted, thus each paramagnetic layer in our 
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structure behaves as a potential barrier for spin-up electrons and a quantum well for spin- 
down ones. Consequently, the spin-up electrons see a double-barrier potential while the 
spin-down ones see a double-well potential. The barrier potential lowers the current density 
of spin-up electrons. For these electrons, the central ZnSe layer behaves as a quantum well, 
thus with increasing the thickness of the central layer, the position of the resonant states 
formed in the well varies and shifts to the lower energy side. This leads to the oscillations 
of the current density which decrease exponentially at each applied bias. Our studies have 
also shown that, with increasing the applied voltage, the transmission coefficients for the 
energies which are near Ep, increase. Since in tunneling process at low temperatures, the 
electrons with energy near Ep carry most of the current, one can see that the current density 
increases with the applied bias. On the other hand, as it is clear in Fig. 1(b), the current 
densities and the amplitude of oscillations for spin-down electrons are much higher than for 
spin-up ones. This is because the structure behaves as a double-well potential, hence the 
transmission coefficients are large even at low energy region and the resonant peaks become 
less sharp. It is important to note that in the spin-down case the resonance is originated 
from the above-well virtual states. 

In Figs. 2(a) and 2(b) we have shown the thickness dependence of the current densities 
for spin-up and spin-down electrons respectively, in an asymmetric (x = 0.04, y = 0.05) 
structure. Here, for x = 0.04, N o:l = 0.27 eV, Vql = —3 meV, T 0L = 1.4 K, and for 
y = 0.05, N a R = 0.26 eV, V 0R = meV, T 0R = 1.7 K Q- At low applied bias, the 
current density and the amplitude of oscillations for spin-up electrons slightly increase when 
the thickness of L2 increases, while with increasing the applied bias a reverse behavior 
appears. An important point to note is that the values of current density for the spin-up 
electrons in the asymmetric structure is three orders of magnitude higher than the symmetric 
one, while in the spin-down case, the values do not change considerably relative to those 
obtained in the symmetric structure. However, the behavior of current density for the two 
spin subbands is completely different in both structures as shown in Figs. 1 and 2. In the 
asymmetric structure and in the spin-down case, with increasing the applied bias and the 
thickness of the central ZnSe layer, the current density increases, which is important for 
possible technological applications. This is a direct consequence of the effect of conduction 
band offset and the difference in the magnetic impurity concentrations in the paramagnetic 
layers. 
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In order to further see the effects of the resonant states in our system, we have shown 
in Figs. 3 and 4 the degree of spin polarization as a function of L 2 under several applied 
voltages in the symmetric and the asymmetric structures, respectively. One can easily 
see in Fig. 3 that, the output current of the symmetric structure exhibits high values of 
spin polarization, which can reach 100% when L 2 increases. Therefore, the results display 
obvious behavior of spin filtering effect in the symmetric structure, which is independent of 
the applied voltage. The spin polarization in the asymmetric structure has been shown in 
Fig. 4. The results show that (in contrast with Fig. 3) in the asymmetric structure the spin 
polarization strongly depends on the applied voltage. In such structure at fixed L 2 the high 
(low) values of spin polarization correspond to low (high) applied voltages. As the thickness 
L 2 increases, the degree of spin polarization for the low voltages decreases slowly, while the 
amplitude of oscillations increases. However, the system for the high voltages exhibits a 
reverse behavior. The results once again indicate that the spin-dependent resonant states 
which are responsible for the oscillations in spin currents, strongly depend on the external 
electric field and the geometrical structure. 

Therefore, the obtained results in the asymmetric structures clearly illustrate that the 
current density and hence the degree of spin polarization can be tuned by changing the 
applied voltage and/or the thickness of the central ZnSe layer. 



IV. SUMMARY 



Based on the transfer matrix method and the effective-mass approxima- 
tion, we have investigated the effect of resonant states on spin transport in 
ZnSe/Zni_ a; Mn a ;Se/ZnSe/Zni_ y Mnj / Se/ZnSe heterostructures. The numerical results 
show that the symmetric structure enables the generation of spin-polarized injection 
currents, since the heterostructure for each value of the applied bias and the thickness of 
the central ZnSe layer, filters out the contribution of spin-up electrons in the total current 
density. On the other hand, our numerical calculations indicate that the asymmetric 
structure can be a good spin filter, if we adjust the applied voltage and the thickness. Such 
behaviors are originated from the enhancement and suppression in the spin-dependent 
resonant states. These interesting features are relevant for devising tunable spin- dependent 



electronic devices such as spin switches, tunable lasers |2J 



241 ]. magnetic resonant tunneling 
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diodes |lOj], and also useful to study fundamental effects involved in the spintronic field. 
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FIG. 1: (a) Spin- up and (b) spin-down current densities as a function of the thickness L2 
ZnSe/Zni-zMn^Se/ZnSe/Zni-zMn^Se/ZnSe heterostructures. 
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FIG. 2: (a) Spin-up and (b) spin-down current densities as a function of the thickness L2 
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FIG. 4: Degree of spin polarization as a function of the thickness L2 for 
ZnSe/Zni_ x Mn a; Se/ZnSe/Zni„j / Mnj / Se/ZnSe heterostructures. 
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